Local polariton states in impure ionic crystals 



V.S. Podolsky, Lev I. Deych and A. A. Lisyansky 
Department of Physics, Queens College, CUNY, Flushing, NY 11367 

(February 1, 2008) 

Abstract 

We consider the dynamics of an ionic crystal with a single impurity in the 
vicinity of the polariton resonance. We show that if the polariton spectrum 
of the host crystal allows for a gap between polariton branches, the defect 
gives rise to a novel kind of local states with frequencies within the gap. 
Despite the atomic size of the impurity we find that new local states are 
predominated by long-wavelength polaritons. The properties of these states 
are shown to be different from the properties of the well-known vibrational 
local states. The difference is due to the singular behavior of the density of 
states of polaritons near the low-frequency boundary of the polariton gap. 
Assuming cubic simmetry of the defect site we consider a complete set of the 
local states arising near the bottom of the polariton gap. 
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I. INTRODUCTION 



An existence of local excitations caused by crystal defects is a well-known phenomenon in 
solid state physics . A single point-like defect can give rise to new states localized in the 
vicinity of a defect with frequencies outside the bands of extended states of a host crystal. 
These states and their interaction with electromagnetic waves (IR absorption, Raman scat- 
tering) were extensively studied since the early forties following pioneering works by Lifshitz 
0-^. In all these studies electromagnetic field was considered as an external field, which 
excites vibrational states. The feedback effect of the vibrations upon electromagnetic field 
was neglected. However, in ionic crystals in the region of frequencies closed to the crossing 
point of phonon and photon dispersion curves one has to take this effect into account because 
the coupling between electromagnetic waves and vibrations becomes so strong that the new 
kind of excitations - polaritons - emerges. In this situation the analyses of local vibrations 
and their interaction with electromagnetic field has to be done more carefully with taking 
the polariton effects into account. For the first time such a consideration was carried out 
in Ref. where a novel kind of local excitations was discovered. It was shown in Ref. ^ 
that if the polariton spectrum of a crystal exhibits a gap between polariton branches, where 
the density of states of a host material is equal to zero, a defect embedded in the crystal 
gives rise to local states with frequencies within the gap. These states are a mix of elec- 
tromagnetic field and excitations of a host material. Their properties appeared to be quite 
different from those of known local phonon states. Unlike pure phonon systems, in the case 
of isotropic materials there is no minimum critical value of the "strength" of the impurity 
for local polariton states to appear. This feature was shown to be caused by a negative 
dispersion of optical phonons resulting in a non-monotonic dispersion of polaritons. 

A general analysis made in Ref. refers to arbitrary "polarization waves." It is applica- 
ble to a variety of excitations, such as phonon-polaritons, exciton-polaritons, etc. Concerning 
a particular case of phonon-polaritons the model of a dipole interaction between scalar elec- 
tromagnetic and polarization waves, used in Ref. |]^, needs to be extended to account for a 
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vector nature of the excitations. Moreover, the analysis of non-isotope impurities, affecting 
the elastic bonds around a defect, is not consistent within the scalar model. 

In the present paper we take into account the vector nature of electromagnetic waves 
interacting with vibrations in crystals with cubic symmetry. We assumed that an impu- 
rity atom in addition to having a different mass can locally change elastic constants. The 
anisotropy of the crystal is assumed to be weak and is neglected in the long wavelength 
limit. We obtained two series of local states which differ in parity. In agreement with results 
of Ref. 0, all these states appear first at the bottom of the polariton gap for infinitesimally 
small variations of impurity parameters. This is shown to be caused by a singularity of the 
density of states in the lower polariton band. This singularity also provides that the local- 
ization of transverse polaritons is most effective near the lower boundary of the gap. We 
show that the local polariton states, unlike usual transverse extended states are affected by 
the interaction with the longitudinal phonon modes. This interaction narrows the frequency 
range available for the local polariton states. 



II. POLARITONS IN A PURE CRYSTAL: ISOTROPIC APPROXIMATION 



We consider below a body centered cubic (BCC) dielectric crystal with two oppositely 
charged ions per each elementary cell. The interaction between ions is assumed to be a 
central one. Denoting masses of the positive and negative ions as m+ and m_, respectively, 
and their displacements as U±(r), where r belongs to the positive or negative sublattice, we 
can write down the equations of stationary lattice vibrations coupled to a coherent electric 
field, Efr): 



m±cj2U±(r) ±gE(r) 



E 



(1) 



where q denotes an ion charge, u is the frequency. 

The right-hand side of Eq. (1) presents all elastic forces acting on the ion hosted at the 
site r. Vectors and denote radius- vectors of neighboring ions in the s-shell spheres of 
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the original and alternative sublattice, respectively. In the case of a central interaction the 
elastic forces within one sublattice have a form: 



f(rR.) 



R,{r,- [U±(r)-U±(r+R3)]} 



(2) 



and forces between ions from different sublattices are: 



f (rR',) = -^^R'. {R'. • [U±(r) - U^(r+R',)]} , 



R 



(20 



where elastic constants of intra- and inter-sublattice interaction, /3{R) and P'{R'), depend 
on a distance between ions only. 

A coherent electric field induced by the ionic vibrations, E(r), invokes an additional pair 
of equations: 



V-E(r)= -47rg^ [u+(l) • V5(r-l)-U_(l+b) • V5(r-l-b) 



(3) 



Vx[VxE(r)] + ^E(r)=-l!^ 
c c 



^[u^(l)5(r-l)-U_(l+b)5(r-l-b) 



(4) 



where vectors 1 denote lattice vectors, b is the basis vector and c is the speed of light. 

The lattice normal modes arise as simultaneous solutions of Eqs. (1,3,4). Fourier trans- 
formation of this system gives the dynamics equations in the momentum representation: 



A_fk) 



= 0. 



(5) 



D(k) + F(k) - uj^m+ D'(k) - F(k) 
\^ D'(k) - F(k) D(k) + F(k) - A. j y 

Here A±(k) denote the Fourier amplitudes of the displacements, D(k) is the dynamical 
matrix of the intra-sublattice interaction: 



(6) 



•-s R' 



and D'(k) is the dynamical matrix of the inter-sublattice interaction: 

D'(k) = - E E R'.®R'.e'''-^'% 

s R' 



(60 



where R(8)R denote the direct products of two vectors. 

The Fourier amphtudes of the field induced by the lattice vibrations can be expressed as 
E = — F (^A_^—A_^/q, with operator F defined as follows: 

where I is the unit tensor, Pji*^ = kakis/k"^ and P"'^ = dap — kakp/k^ are the longitudinal 
and transverse projectors, respectively. 

The dynamical matrices given by Eqs. (6,6') are of a general type and solutions of 
Eq. (5), in general, do not split into longitudinal and transverse modes. However, consid- 
ering long-wave excitations in weakly anisotropic crystals we can make use of the isotropic 
approximation for the dynamical matrices: 

D(k)=7||(fc)P||+7x(fc)P±, (8) 

D'(k)=7[|(fc)P|l+7^(fc)P^, (8') 

where scalar functions, 7o.(/c) and ^^ik), can be expressed in terms of frequencies of lon- 
gitudinal and transverse phonons. In a crystal of cubic symmetry the dynamical matrices 
become trivial at the center of the Brillouin zone: 

D(o)= -mo)=E^j:K^K-\ii:f^'iR's)z's, (9) 

where Z'^ denotes a total number of ions in s-shell and we make use of the identity, 
X^R/ Rg (8) Rs = IZ^i?^^/3, valid in crystals with cubic symmetry. Eq. (9) sets a condition 
for parameters in Eqs. (8), so that all 'j„{0) and —7^(0) are equal to a positive constant 

In the isotopic approximation all normal modes of a crystal become either longitudinal 
or transverse, A^^\k) = e„{k.)A^^\k), where e^Ck) are longitudinal or transverse unit polar- 
ization vectors, cr is a polarization index. Eqs. (5,8) give the following relation between the 
Fourier amplitudes of displacements in different sublattices: 
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^ fa la (10) 

la + f^-uj^m+ 

For longitudinal modes /y =47rg^/a^ = /. The corresponding dispersion equation, 

-fc^ = u-u^i..) u-u.^ „) = 0, (11) 

defines acoustic, c<jy(/i;), and optical, fiy(/c), branches of longitudinal phonons. For transverse 
modes f^ = fuj'^/ (u'^ — c^k^) , therefore, the dispersion equation, 

y m+ J y m_ J m+m^ tu'^ — c^k'^ ' 

gives one acoustic, u;^(fc), and two polariton, f2^(fc), branches. 

Neglecting the field effects and solving Eqs. (11,12) in the long wavelength limit one can 
obtain expressions for '~fa{k) and 7^(A;) in terms of conventional parameters: 

'y,{k)^-f+fie{vl-v'^'), (13) 



(13') 



'M 

where n and M are reduced and total masses of ions within the elementary cell, respectively, 
V(j and f ^ are the velocities of acoustic and optical phonons with a given polarization. 

Eqs. (11,12) show that the internal field affects dispersion relations of all lattice excita- 
tions. However, the physical effects of the photon-phonon interaction are substantial in a 
vicinity of the polariton resonance which takes place in the long wavelength region. It is 
straightforward to show that the acoustic branches, 

ulik) ^ 2^^ ^ vie, (14) 

are unaffected by the field, whereas, the interaction with the field results in the uniform 
up-shift of the longitudinal optical branch, 

uj^=nl{k)+d^, (15) 

and leads to the well-known polariton dispersion relation for optical transverse excitations 
ii, 



= dV. (16) 
Eq. (16) describes two polariton branches, D,±{k), with corresponding dispersion laws: 

nl{k)^ 1/4 (^'J[n4k)+ckf+d'^±^ [n^k) -ckf+ ^ , (i7) 

where D,± and ily are frequencies of the transverse and longitudinal optical phonons, respec- 
tively, 

nUk)^^-2^1^ ^u:l-v^'k\ (18) 

a phonon- photon coupling parameter, (P = f / fi = inq"^ / fj^a^ , is the ionic "plasma frequency" 
and u!q = j/ij, is the optical activation frequency 

Analysis of Eqs. (14,15,17) shows that the lower polariton branch, Q^(A;), has a zero 
activation frequency and extends over both acoustic bands, whereas, the longitudinal optical 
branch, fl'^^{k) + (P, overlaps with the top part of the polariton gap. Therefore, a truly 
forbidden gap, with no modes of any kind inside, may exist only between the lower boundary 
of the polariton gap and the bottom of the longitudinal band [Fig. 1]. When the dispersion 
of optical phonons is neglected, both polariton branches are monotonic and the spectral gap 
is between the transverse and longitudinal frequencies, cUq and cUq + (P. 

It turns out that accounting for the phonon dispersion causes a qualitative change of this 
picture. First, the upper boundary of the frequency gap, a;|, is now set by the bottom of 
the optical longitudinal band. Second, the lower polariton branch in the case of a negative 
dispersion becomes non-monotonic and gains a maximum at some k = ko, close to the 
center of the Brillouin zone [Fig 1.]. Calculations in the long wavelength approximation give 
/cq ~ 2u!od/v'j_c and a new lower boundary of the gap, cuf — ^}'i,{ko) ~ a;Q— 2v'j_u!od/c. Since 
2v'j_uJod/c ~ v'lk"^ <^ lOq, the result is consistent with the approximation. However, because 
c^/cq ~ 2cu!od/v'j_ ^ lUq, the maximum of fl'^{k) is far away from the very narrow (due to 
c ^ v'j_) polariton resonance region. Since Q^(A;) reaches its maximum at the surface of a 
finite area inside the Brillouin zone, the density of polariton states diverges at the gap's lower 



UJ 



-nl{k) 



CO 
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boundary. As we show latter, it causes an absence of a lower threshold for local polariton 
states. 



III. POLARITON LOCAL STATES 



When a host ion at the site r = of the positive sublattice is replaced by an impurity ion 
with the same charge, it causes a local deviation of the crystal density and a local change of 
elastic constants. To account for these facts we need to add extra forces to the right hand 
side of Eq. (1). For positive and negative sublattices, respectively, these forces are: 

5f+(r) = |cu25mU+(0) - ^ n [n ■ U+(0) -n ■ U_(n)] | 5ro (19) 

and 

5f-(r) = -^^n[n-U+(0)-n-U_(n)l5r„, (19') 

where Sm is the difference between the masses of an impurity and a host ion, 5(3 is a shift in 
the elastic constant in the impurity's near-neighbor shell and vectors n denote radius-vectors 
of the impurity's nearest neighbors. 

The dynamic equation (5) is now modified: 



'^D(k) +F(k) -cu2m^ 



D'(k) - F(k) 
\^ D'(k)-F(k) D(k) +F(k) 
where N is the number of ions in one sublattice and 



^^A.(k)^ 



A_(k) 



1 

N 



(20) 



B_ 



-zk-n 



(21) 
(21') 



Here U_,_ = U_,_(0) and all U_(n) appear only in combinations C/„ = n-U_(n), because 
ion-ion forces are central. 

Let us denote the matrix in the left-hand side of Eq. (20) as Cartesian {a) and sublattice 
(e) indices of the amplitudes A'^. In the isotropic approximation this matrix can be decom- 
posed, L = L|| P|| + ® Pj^, where 2 x 2-matrices L^, operating on sublattice indexes 
only, are defined as the follows: 
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''a ''a 



The inverse matrix, L ^ = Ly ^ 



\ ia- i<y la + fa- UJ^m_ J 



(22) 



-1 



^ 9t g'a ^ 



G = G|| (g) P|| + G_L ® P^, where 



la + fa {la-f<^) 



771- 



m+rri- 



y-l'a It J 



(23) 



\9a 9a J 

is the Green's function of the system in the momentum representation. In accordance with 
a general property of Green's functions, poles of its matrix elements, ga{uJ, k), coincide with 
eigenfrequencies of a pure crystal. 



g^{uj) oc (uj^-ujI) ^io'^-nl) \uj^-nl) \ 

as it follows from Eqs. (11,12). 

Solving Eq. (20) with the help of Eqs. (22,23), one can obtain Fourier amplitudes of the 
ion displacements, A+(k) and A_(k), and the displacements themselves: 



U+(r) = E / e*'^-e<,(k)®e<,(k) (^+B+ + ^;B_, 
U_(r) = E / e*'^"e.(k)®e.(k) UB^ + g;B_, 



(24) 



(240 



where a symbol /(j^.) denotes (l/A^)I]k ~ ('^/27r)^/ dk, with wave vectors taken from the 
first Brillouin zone only. 

Substitution of B+ and B_, given by Eqs. (21,21'), into Eqs. (24) allows one to express 
all displacements in terms of U_,_ and only, and then obtain a closed system of equations 
for these variables. In the case of an isotope impurity, considered earlier in the scalar model 
of Ref. 0, we set 5/3 = and this system of spectral equations reads: 



6mJ2 [ ^+e,(k)®e^(k)U+. 

a "'(k) 



(25) 
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In a cubic crystal, vector U+ can be arbitrary [Appendix A], whereas frequencies of excita- 
tions are defined by the equation: 

1 = ^ (4 + ^91) = Smco'lico^ (26) 

which generahzes Eq. (10) obtained in Ref. [0. 

All solutions of these equations can be divided into two classes: extended and local 
states. Since an impurity destroys the translational symmetry of the crystal, any state is 
now a superposition of all normal modes available in the first Brillouin zone. For extended 
states, corresponding to scattering states with well-defined wave vectors, frequencies fall into 
the bands of a pure crystal. Local states, dependent upon the value of the parameter 6m 
in Eq. (26), may arise outside of the bands. From the structure of g^lujk.) it might seem 
that when a frequency is close to a particular band, the modes from it dominate in the 
corresponding state. However, because contributions of the near and distant bands could be 
weakened or strengthened by the low or high density of states in them, direct calculations 
are required here. 

As we already mentioned, the function Q'^{k) reaches its maximum value, ul, at k = ko, 
close to the center of the Brillouin zone. Its expansion around ko does not contain a linear 
term: 

Q^_{k)=ujl-u^a^ {k-kof + O [a^ [k-k^f] . (27) 

It immediately shows that the integral in the left-hand side of Eq. (26) diverges at uj = uji 
due to a contribution from the lower polaritons branch. Therefore, considering uj in the 
frequency gap close to its lower boundary we can omit g'^ from Eq. (26) and rewrite J^y^ g~\_ 
using the density of states, P-(e), in the lower polaritons band: 



where -F(fc) denotes all factors of g\^ other than (^cu^ — fi^^ and the last integration is 
performed on a complex e-plane [Appendix B]. 
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Near the bottom of the polariton band, for small e, the density of states has the usual 
Kohn shape, p^{s) oc ^/6. The asymptote of P-{e) near the top of the band, for small 
z = ujf—e, can be found with the help of Eq. (27): 

P-^ ' \2ttJ J \Vnl{k)\^ \2ttJ 27r2z/2V2' ^^^^ 

where all omitted terms are regular at z = 0. 

Since for u^uJi + O the divergent part of is provided by a region of small z, in 

Eq. (28) we can replace the exact density of states with the found asymptote. Considering 
u'^—uf <^ ujf it allows us to calculate the leading part of J+(co'^) and transform the spectral 
equation (26) into the form [Appendix B]: 



m+ 2 {akQ)"^ uJq m_ 



Using estimates made in section. 2 and recalling that ai/ = f^, we obtain: 



(30) 



[31) 



m a^ujQd 



This result recovers Eq. (14) of Ref. obtained in the scalar model and supports the 
conclusion of the absence of the lower localization threshold for light isotope impurities. 

In a general case of non-isotope impurity the spectral system is of the 11th rank with 
11 variables, U+ = U+(0), f/„ = n • U_(n). The system can be further simplified with the 
help of the crystal symmetry. The exact point group of a cubic crystal includes the space 
inversion. Therefore, all excitations can be classified by their spatial parity. 

For the odd states, where U^(— r) = — U_|_(r), one can see that both U_,_ and B_,_ are equal 
to zero and the rank of the spectral system reduces from 11 to 4. Considering displacements 
of four negative ions at corners of one face of the BCC-lattice elementary cell (see Fig. 2) as 
independent variables, from Eqs. (24) we obtain the following system: 

4 2 

J2M,,,U,= - — Us, (32) 
s'=i 25/5 

where index s numerates the chosen ions, their radius-vectors are denoted as n,, and Ug = 
=ns • U_(ns). The matrix Mgs' here is given by the equation: 
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Mss' = M{ns, ns>) = ^ 9a{^T- nj(e^ • n^) sin(k • n^) sin(k • n^/). (33) 

Using transformation properties of the integrand one can show that the matrix Mgs' remains 
invariant, M(Qn, Qn') = M(n, n'), under any point transformation Q. Therefore, the 
elements of the symmetric matrix Mss' are equal for equivalent pairs of ions: Mu — M22 — 
M3S — M44, then M12 — M23 — M34 and Mi3 = M24. After finding eigenvalues of this matrix 
we obtain three spectral equations: 

/xo = Mn-Mi3, (34) 



25(5 

A(± = Mii + Mn±2Mi2. (34') 



25/3 

The structure of the eigenvectors shows that the states corresponding to /i+ and represent 
"rhombic" and "tetragonal" oscillations localized around the stationary impurity, whereas, 
the states corresponding to involve both types of deformations of the elementary cell 
[Appendix C]. 

Near the lower boundary of the gap the arguments used in evaluation of Eq. (26) are 
also applicable. Retaining the leading terms in Eq. (33) and using symmetry properties of 
the arising expressions [Appendix C], we obtain the following expressions: 

1 

o 



n^-(ni-n3)^ / k'gl, (35) 



if 

At±~ -Uni-n3)2 + y ±2 



4 

(ni -na) - — 



From the geometry of the elementary cell it follows that (rii • n2) = — (rii ■ n3) =72^/3 = 0^/4 
and, therefore, /Iq — /i^ — (o^/lO) k^g±- It leads to the unique spectral equation for all 
"tetragonal" modes [Appendix B]: 



6(3 ^ 15nu^u^-uU My 

This result shows that the odd local states arise upon an infinitesimally small strengthening 
of local bonds associated with the impurity, 5/3 > +0. This effect, similar to the isotope 
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impurity case, is due to a singularity of the density of states in the lower polariton band. 
Eq. (36), compared to Eq. (30), has an additional small factor, (a/co)^ • Therefore, for the 
same relative deviations, Sp/ j3 and 5m/ the odd states lie much closer to the gap's bound- 
ary than the states associated with an isotope impurity. Moreover, these states are much 
less sensitive to variations of parameters than the isotope-induced local states. This fact 
reveals itself in an extreme form for "rhombic" modes corresponding to Eq. (35') gives 
//+ = which means that within the approximations used these states appear right at the 
gap's bottom for any value of 513. Accounting for the higher-order terms in expansions of 
the sin -functions in Eq. (31) will separate these modes from the gap's boundary but the 



separation, ^Ju^—u^ oc 5(5 {akof , remains the smallest among all considered local states. 

For the even states, where U_|_(— r) = U_|_(r), the spectral system contains 7 independent 
variables. Those are components of the impurity displacement, U_,_, and radial projections 
of the displacements of four chosen neighboring ions, Ug- The elements of the matrix in the 
corresponding spectral system can be written as integrals of the Green's function elements, 
similar to Eq. (33). Considering local states near the gap's bottom, we again can retain in 
these integrals only transverse terms and use the power expansions of non-singular factors 
of the integrands. As it was shown above, /(k)(o^)^9'I <^ I(k)9± due to the dominant 
contribution of the small-A; region in these integrals. Therefore, since the local states near 
the gap's bottom are, roughly, "made" of long wavelength polaritons, we can disregard the 
exponential factors in the matrix elements of the spectral system [exp(±ik • n) 1] and 
neglect in Eqs. (21) differences between displacements of identical ions within the elementary 
cell, U_ ~ U_(n). These approximations lead to the following spectral system for even 
states: 



U 



\ 



( 



J 



V 



u'^5m-S5(5/2> 85p/3 
85P/3 -85(3/3 



\ 



r r 

The corresponding spectral equation reads: 

85/5 .2c___r+,. 8uj'^5p5m 



U_ 
U. 



(37) 



(/++/- -2/') -cu25m/++ 1 



(38) 
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where all /-factors are straightforward to evaluate near the gap's bottom [Appendix B]: 



^-"^^ (39) 



^ 3 /k) ^"^ ~ Q^„, /,.,2_,.,2 m.m^ ' ^^^"^ 



The right-hand side of Eq. (38) is proportional to a determinant of a degenerate operator, the 
transverse propagator, and it must be equal to zero identically. Substitution of Eqs. (39,39') 
into Eq. (38) shows this explicitly and transforms the spectral equation of the even states 
into the following one: 



2 {akof 



Snufji'^uj^—uji L 



8SP n / m 
— UQdm 



M 



1. (40) 



One can check that here ions displacements satisfy the relationship, m?_U+ -|- m^U_ = 0, 
corresponding to optical vibrations. 

If we set 5/3 = in Eq. (40) it reproduces Eq. (30) obtained earher for an isotope impurity. 
That equation has a solution only when 6m < 0. In a general case even local states exist 
only if variations of parameters satisfy the inequality: 

5P-lu;l5m(^^)\o. (41) 

Eqs. (36,40) allow us to outhne regions of the local states [Fig. 3] on a plane of im- 
purity parameters, 5m, 5/3. Taking into account obvious physical limitations, 5m > — m+, 
5P > — /3, one can see that odd states appear in the right upper quadrant bounded by the 
lines 5m — —m+, 5(3 — 0. The region of even states is to the right of 5m — —m+ and above 
the critical line, 5/3 — (3/8) {m- / M)^ UQ5m. All other areas of {5m, 5/3)-plane are blocked for 
the polariton localization. 

Since we deal with local states near the gap's bottom, our analysis is valid in a vicinity 
of the critical hne for even states, and close to 5m-axis for odd states. When the impurity 
parameters move outside these regions the frequencies of the corresponding local states move 
away from the gap's bottom. The lines where the frequencies approach the top of the gap 
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establish the outer boundaries of the locahzation regions. Unhke the situation near the 
gap's bottom, all terms of all spectral equations remain finite when uj tends to uj2i what 
guarantees an existence of limited localization regions. On the other hand, this makes it 
impossible to do any rigorous analytical calculations. However, using the density of states 
representation of the spectral equations, analogous to Eq. (28), and utilizing some trial 
functions to simulate densities of states in all bands, one can obtain qualitatively reliable 
results. Our preliminary estimates show that near the upper boundary of the gap local states 
are substantially composed of transverse and longitudinal phonons. The balance between 
them depends on widths of the corresponding bands and the polariton gap. The contribution 
of the long wavelength transverse polaritons into these states is proportional to (f /c)^ and is 
negligible. Roughly, it is caused by the fact that the density of states is inversely proportional 
to for phonons and to cf ^ for long wavelength polaritons. A more detailed analysis of 
states located far away from the gap's bottom will be presented elsewhere. 



IV. DISCUSSION AND CONCLUSIONS 

We have considered local polariton states in BCC ionic crystals. It was assumed that 
the crystal anisotropy is weak and can be neglected in the long wavelength limit. This 
approximation was proved to be self-consistent for states located near the bottom of the 
polariton gap. We found two series of local states, different in parity. The new states 
appear right at the bottom of the polariton gap upon infinitesimally small variations of an 
impurity parameters. This is in contrast with 3-D phonon systems where a lower threshold 
for local states always exists. In Ref. |^ the general theorem regarding the presence of 

the threshold for arising local states in bandgaps of periodic systems were given. However, 
the proof assumed the finite values of density of states in the entire band of pure system. 
We show that the singularity in the density of states in the lower polariton band causes 
the absence of this threshold. This singularity also provides that the states near the gap's 
bottom are formed mostly by the long wavelength transverse polaritons. The local states 
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move toward the upper boundary of the gap upon increase of impurity parameters Sm and/or 

We have outhned regions of new local polariton states on a plane of impurity parame- 
ters. In a region where both types of states coexist, the odd states precede the even ones. 
They appear first at the gap's bottom and remain near it when the impurity parameters 
vary much longer than the even states. Local states are a mixture of transverse phonons 
and photons. Comparing amplitudes of the field, E, and the crystal polarization, P, one 
can make an estimate of the energy partitions of local polariton states. Because the char- 
acteristic momentum, ko, happened to be away from the polariton resonance region, the 
ratio E/P turns out to be of the order of v'j_/c. An account for the electronic polarization 
of ions renormalizes c and increases this ratio. More detailed analysis of these aspects of 
the polariton local states will be done elsewhere. Condition (41) can help in a search for 
compounds where local polariton states can be observed. 

The results of this paper were obtained within the harmonic approximation of crystal 
dynamics. Phonon-phonon interaction caused by a lattice anharmonism will affect this 
picture. Phonon-phonon scattering causes damping of all elementary excitations and the 
local states as well. It leads to a broadening of all spectral lines and it also washes out 
all sharp boundaries in the initial excitation spectra. However, it seems physically evident 
that the scale of this broadening is far below the typical width of phonon bands and cannot 
change the topology of the initial band structure, or close the polariton gap. If we dress 
up the elementary excitations and renormaUze their spectra, the maximum of the density 
of states in the lower polariton band and the corresponding singularity will persist after 
the renormalization since it has a topological origin. That, in turn, guarantees the absence 
of a threshold for local states near the (renormalized) gap's bottom at zero temperature. 
However, at non-zero temperature thermal broadening of spectral lines will set a finite 
threshold for local states. 

Another factor that leads to a threshold is crystal anisotropy. It provides a difference 
between maxima of (k) in different directions that causes the density of states in the lower 
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polariton branch to be finite everywhere. It should be emphasized, however, that properties 
of the local polariton states remain quite different from the corresponding properties of 
purely phonon local states. Particularly, one can see from our results, that eigenfrequencies 
of local polaritons are much more sensible to the crystal structure of the host materials. 
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It is known that under any point transformation, Q, the frequency of any normal mode 
remains invariant, a;^(Qk) =a;^(k), and the corresponding polarization vector transforms as 
the follows: e(Qk) = Q~-'^e(k). Because the Brillouin zone also maps exactly into itself, one 
can see that 

t= / g[uj\k)]e(k)®e(k)= [ ^/b^fQk)] e(Qk)®e(Qk) = Q-^tQ. (42) 

>/(k) ^(Qk) 

Therefore, any tensor of this type is an invariant of the point group of a crystal. In the 
cubic system any group-invariant operator must be trivial, T — tl, since the group contains 
non-coUinear axis of different order. Calculating the trace of the operator, one can find: 

t = l/3 / o[u;2(k)][e(k)-e(k)] i, (43) 
J{k) 

where I is the identity operator. 

The density of states in the lower polariton band, P-{£), is defined at the complex e- 
plane, cut from to o;^. The contour of integration, C, in Eq. (26) runs along the upper side 
of the cut (it corresponds to integration over k < ko), turns around its right edge and then, 
for the integration over k > kg, returns along the lower side of the cut to some point, uj'i, 
fixed by the bottom of the optical photon band. 

As it was discussed in the text, the leading part of /(ciJ^) at A(tc;) =tu^— cj^ ^ ujj, comes 
from the region close to the right edge of the cut, e—u;f. Therefore, evaluating Eq. (26), we 
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can use there the asymptote of the density of states given by Eq. (27). In addition, as it 
follows from Eq. (25), k—ko— ^z^^'^/au near the right edge of the cut at its upper and lower 
sides, respectively. Taking this into account when expanding F{k) around ko in Eq. (26), 
we obtain: 



+ . 2[Fo-F^,/^/aiy+---]dz r'-lW 2 [Fo + F^,/^/aiy+ ■ ■ ■] dz 

' Jo 3^{A + z) Jo 3Vi(A + z) 



(44) 



A rescaling z zA shows that contributions in I~^{uj^) from all terms of expansion of 
F{k) tend to zero in the hmit A— >+0, except for the first two terms. Further elementary 
integration gives the singular part of the considered integral: 

2{ako)^Fiko) _ 2{ak^f {uj^ -^kl)\^^{k^) + f^{k^)-m.u^] 



(45) 



37ri/AV2 37ri/AV2 m+m-[uj^-nl{ko)][uj^-ujl{kQ)] ' 
This result is asymptotically exact since all omitted terms are regular at A — > +0. Using 
estimations made in Section 2. and taking into account conditions, ako <^ 1, v'^kg <^ Uq, 
(P'kQ ^ uJq and a; ~ a;i ~ ujq, it is straightforward to obtain: 

2{ak^f 



UJ^ — UJ 



2 Mm 



(46) 



+ 



which immediately leads to Eq. (28). 

In a similar way one can obtain Eqs. (37), as well as, calculate the integral in Eqs. (33): 

2 {ako)'^ uJq m+ 



2 - ^ 2(afeo)^ {u;'-cX) bAko) + f±{ko)-m+u'] ^ 



(47) 



Eq. (30) with the matrix Mgs' defined by Eq. (31) can be rewritten as follows: 







































Mi3 




Mn+Mi3-/i 


2Mi2 
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Mi2 


Mi3 


2Mi2 


Mn + Mi3 







(48) 



/ 

which has three eigenvalues given by Eqs. (32). 

For fi = fi± this equation gives Ui = U2 = and U3 = ±1/4. Recalling that Us = 
n • 11(115) and the states we consider here are antisymmetric, one can see that in /x+-mode 
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the elementary cell inclosing a defect experiences "rhombic" deformations. In //_-mode, 
since the ion motion is anti-phased, it produces "tetragonal" deformations. For fj, = fj,Q 
Eq. (CI) leaves Ui and U2 independent, whereas, U3 = Ui/2 and t/4 = t/2/2. It results in 
combined, "tetra-rhombic" deformations of the elementary cell in these modes. 

To obtain the explicit expressions of the eigenvalues, /x, we have to calculate the matrix 
elements Mgs' ■ When the frequency of the considered states is close to uji we can retain only 
first terms in expansions of sin -factors in Eq. (29). This approximation is self-consistent, 
because the region of small wave vectors is the major contribution in Eq. (29). Carrying out 
the calculations, one needs to know tensors of two types: f{k)k"'k^ and f{k)k"k'^k"k'^, 
where f{k) is an arbitrary invariant function. Prom Appendix A it immediately follows 
that, f^f{k)k°'k^ — {1/3)6"^ J^f{k)k'^. The second tensor is obviously totally symmetric, 
therefore: 



where summation over a = P and a = (3 defines the numerical factor, / = 1/15. Using 
these results in Mg^/ and retaining there the transverse terms only, we obtain: 



/ f{k)k'^k^k''kf^ = I (^"'^r'^ + S^^S^^f^ + S^^^S^") j f{k)k\ 




(49) 



ss 



1 




+ 2(n,n,0^ 
5 



) 



I. 




(50) 



that, in turn, leads to Eqs. (33). 
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Figure Captions 



Fig. 1. Elementary cell of a body centered cubic lattice. Here a is a lattice parameter, 
indexes s, s' numerate ions within the near-neighbor shell, iris are radius- vectors of the ions. 

Fig. 2. Phonon-polariton dispersion curves in the isotropic approximation. Here lines 
cu — D,-i-{k) present the lower and upper polariton branches, D,±{k) and u!+\\{k) are transverse 
and longitudinal optical phonon branches, uj±{k) and a;_||(A;) are the transverse and longitu- 
dinal acoustic phonon branches, respectively. The horisontal lines cuq and cus are the upper 
and lower boundaries of the transverse optical phonon band, ^Jujq-\- d? and UJ2 are the upper 
and lower boundaries of the longitudinal optical phonon band, respectively. The polariton 
gap is bounded by lines uji and uj2- The wave number /cq corresponds to the maximum of 
the lower polariton branch. 

The figure presents qualitative features of phonon-polariton spectra in the long wave- 
length region and does not reflect an actual scale of wave velocities. 

Fig. 3. Regions of existence of local polariton states. 5m is the difference between 
masses of the impurity and a host atom, (5/3 is the impurity induced deviation of the elastic 
constant within the near-neighbor shell, m+ is the mass of a host ion in the "positive" 
sublattice, (5 is the elastic constant of a near-neighbor interaction in a perfect crystal. Odd 
local polariton states exist in the region bound by the sohd lines 5m = — m+, 5(3 = 0. The 
region of even local states is to the right of the solid line 5m = — m+ and above the dashed 
line, 5/3 = (3/8) {m^/Mfujl5m. 
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